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Abstract. We present new calculations of the attenuation of stellar light from spiral galaxies using geometries 
for stars and dust which can reproduce the entire spectral energy distribution from the ultraviolet (UV) to the 
Far- infrared (FIR)/submillimeter (submm) and can also account for the surface brightness distribution in both 
the optical/Near-infrared (NIR) and FIR/submm. The calculations are based on the model of Popescu et al. 
(2000), which incorporates a dustless stellar bulge, a disk of old stars with associated diffuse dust, a thin disk 
of young stars with associated diffuse dust, and a clumpy dust component associated with star-forming regions 
in the thin disk. The attenuations, which incorporate the effects of multiple anisotropic scattering, are derived 
separately for each stellar component, and presented in the form of easily accessible polynomial fits as a function 
of inclination, for a grid in optical depth and wavelength. The wavelength range considered is between 912 A and 
2.2 fim, sampled such that attenuation can be conveniently calculated both for the standard optical bands and for 
the bands covered by GALEX. The attenuation characteristics of the individual stellar components show marked 
differences between each other. A general formula is given for the calculation of composite attenuation, valid for 
any combination of the bulge-to-disk ratio and amount of dumpiness. As an example, we show how the optical 
depth derived from the variation of attenuation with inclination depends on the bulge-to-disk ratio. Finally, a 
recipe is given for a self-consistent determination of the optical depth from the Ha/H/3 line ratio. 

Key words. Galaxies: spiral - (ISM) dust, extinction - Radiative transfer - Galaxies: ISM - (ISM) HII regions - 
Galaxies: bulges 



1. Introduction 

The measurement of star-formation rates and star- 
formation histories of galaxies - and indeed of the universe 
as a whole - requires a quantitative understanding of the 
effect of dust in attenuating the light from different stel- 
lar populations. Traditionally, the effect of dust has been 
quantified by statistical analysis of the variation of op- 
tical surface brightness with inclination (see Byun et al. 
1994 and Calzetti 2001 for reviews). These studies often 
reached conflicting conclusions about the optical thick- 
ness of galactic disks, and the need to solve this puzzle 
prompted the development of models for the propagation 
of light in spiral disks using radiation transfer calculations. 

Send offprint requests to: Richard.Tuffs@mpi-hd.mpg.de 



The first work of this type was that of Kylafis & 
Bahcall (1987), who modelled the large scale distribu- 
tion of stellar emissivity and dust of the edge-on galaxy 
NGC 891 using a finite disk and incorporating anisotropic 
and multiple scattering. Using the same method, Byun et 
al. (1994) produced simulations of disk galaxies of vari- 
ous morphologies and optical thicknesses, while Xilouris 
et al. (1997, 1998, 1999) apphed the technique to fit the 
optical/Near-infrared (NIR) appearance of further edge- 
on galaxies. Another approach was taken by Witt et al. 
(1992), who studied the transfer of radiation within a vari- 
ety of spherical geometries, also incorporating anisotropic 
and multiple scattering. This work was extended by Witt 
& Gordon (1996, 2000) to include random distributions 
of dusty clumps within or around a smooth distribution 
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of stars and further extended by Kuchinski et al. (1998) 
to model exponential disks and bulges of observed spiral 
galaxies (see also an application of this technique to ex- 
tremely red galaxies by Pierini et al. 2004). A physical 
approach to dumpiness in terms of its relation to inter- 
stellar turbulence and its effect on attenuation has been 
recently presented by Fischera et al. (2004). 

Specific studies of the attenuation of the spatially in- 
tegrated light from spiral galaxies have been made by 
Bianchi et al. (1996) and by Baes & Dejonghc (2001). Both 
works describe the influence of scattering and geometry on 
the derived attenuation as a function of incUnation, optical 
depth and wavelength. Building on the work of Bianchi et 
al. (1996), Ferrara et al. (1999) presented a grid of atten- 
uation values corresponding to a range of model galaxies 
spanning different bulge-to-disk ratios, opacities, relative 
star/dust geometries and dust type (Milky Way and Small 
Magellanic Cloud type). 

More powerful constraints on the geometrical distribu- 
tions of stars and dust can be obtained by a joint consid- 
eration of the direct starlight, emitted in the ultraviolet 
(UV)/optical/NIR, and of the starlight which is rcradi- 
ated in the Far- infrared (FIR)/submillimeter (submm)^. 
In particular the study of Popescu et al. (2000) showed 
that, in addition to the single exponential diffuse dust 
disk that was used to model the optical/NIR emission, two 
further sources of opacity are needed to account for the 
observed amplitude^ and colour of the FIR/submm emis- 
sion. Firstly, a clumpy and strongly heated distribution of 
dust spatially correlated with the UV-emitting young stel- 
lar population is required to account for the FIR colours. 
This clumpy distribution of dust can most naturally be 
associated with the opaque parent molecular clouds of 
massive stars within the star- forming regions. It is also 
an integral element of the model of Silva et al. (1998) 
for calculating the UV-submm spectral energy distribu- 
tion (SED) of galaxies, and is not to be confused with 
the randomly distributed clumps introduced by Witt & 
Gordon (1996). The second additional source of opacity is 
diffuse dust associated with the spiral arms, which was ap- 
proximated by a second exponential dust disk having the 
same spatial distribution as the UV-emitting stellar disk. 
This component was needed to accoimt for the amplitude 
of the submm emission. Strong support for the reality of 
the elements of the model of Popescu et al. (2000) comes 
not only from the success of the model to fit the inte- 
grated SED, but also from the good agreement between 
the observed surface brightness distributions in the FIR 
and those predicted by the model (Popescu et al. 2004). 
Additional support is provided by the ability to predict 
UV magnitudes, as shown in this paper. 



In the first paper of this series (Popescu et al. 2000; 
hereafter Paper I) we presented the general technique^ 
for modelling the optical/NIR- FIR/submm SED (sec also 
Popescu & Tuffs 2002b) and applied it to the edge-on 
galaxy NGC 891. In the second paper (Misiriotis et al. 
2001; hereafter Paper II) the same model was success- 
fully applied to further edge-on spiral galaxies (NGC 5907, 
NGC 4013, UGC 1082, and UGC 2048), confirming that 
the features of the solutions for NGC 891 arc more gener- 
ally applicable. In the present paper we adopt these fea- 
tures as the basis for calculating the attenuation of stellar 
light, taking into account the constraints on the geome- 
try of stars and dust arising from our fits to the optical- 
submm SEDs. These calculations are made for a grid in 
disk opacity and inclination. We also generalise the model 
so that it is applicable to giant spiral galaxies with dif- 
ferent amounts of dumpiness and with different bulge-to- 
disk ratios (Hubble types). We emphasise that the model, 
in the form presented here, cannot be expected to work 
for small low luminosity systems like dwarf galaxies, which 
have systematically different geometries. Furthermore, the 
composition of dust in dwarf galaxies may differ from that 
in giant spirals, resulting in a different extinction curve 
from the one found for NGC 891 and the other galaxies 
modelled in Paper II. 

A companion paper will present a corresponding grid 
of calculations for the FIR/submm emission, thus pro- 
viding a library of solutions of the SEDs over the entire 
UV-submm range. This paper is organised as follows: In 
Sect. 2 we present tlic model used for the calculation of 
attenuation. Sect. 3 describes the use of radiation transfer 
calculations to obtain the attenuation in the diffuse disk, 
thin disk and bulge components. In Sect. 4 we discuss and 
compare the attenuation characteristics of the disk, thin 
disk and bulge. A formula and a recipe to calculate com- 
posite attenuation of the integrated emission from galaxies 
as a function of the parameters inclination, central face-on 
B-band optical depth, dumpiness and bulge-to-disk ratio 
is given in Sect. 5. In the same section we also discuss 
how to choose values for these parameters and we give an 
example for the calculation of attenuation for the case of 
NGC 891 seen at different viewing angles. To illustrate the 
use of composite attenuation curves in the interpretation 
of optical/NIR data we consider in Sect. 6 the effect of a 
varying bulge-to-disk ratio on the inclination dependence 
of the apparent optical emission from galaxies. In Sect. 7 
we give our summary and conclusions. Readers interested 
only in practical applications of the model are directed to 
Sect. 3 and 5, in particular to Tables 4-6 and Eqs. 6, 17 
and 18. 



^ Based on recent ISO observations of a complete opti- 
cally selected sample of normal late-type galaxies (Tuffs et 
al. 2002a,b) , the percentage of stellar light reradiated by dust 
was found to account for ~ 30% of the bolometric luminosity 
(Popescu & Tuffs 2002a). 



^ A simplified version of this technique has been recently 
applied by Misiriotis et al. (2004) to fit the FIR SEDs of bright 
IRAS galaxies. 
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2. The model 

In this section we describe the characteristics of the model 
adopted for the calculation of attenuation. In terms of 
geometry, the model can be divided into a diffuse compo- 
nent and a clumpy component, the latter being associated 
with the star-forming regions. Direct evidence for the ex- 
istence of these two geometrical components comes from 
recent observations of resolved nearby galaxies done with 
the ISOPHOT instrument (Lemke et al. 1996) on board 
the Infrared Space Observatory (ISO: Kcsslcr et al. 1996). 
For example the maps of M 31 (Haas et al. 1998) and of 
M 33 (Hippelein et al. 2003) clearly show a diffuse disk of 
cold dust emission prominent on the 170 /im images as well 
as warm emission from HII regions along the spiral arms, 
prominent on the 60 /um images. Warm and cold emission 
components have also been inferred from statistical stud- 
ies of local universe galaxies observed with ISOPHOT, in 
particular from mapping observations covering the whole 
optical disk (Stickel et al. 2000 for serendipitously detected 
spirals; Popescu et al. 2002 for Virgo cluster galaxies). 

2.1. The diffuse component 

The diffuse component is comprised of a diffuse old stellar 

population and associated dust and a diffuse young stel- 
lar population and associated dust. The diffuse old stel- 
lar population has both a bulge and a disk component, 
whereas the diffuse young stellar population resides only 
in a thin disk. Throughout this paper we will use the su- 
perscript "disk", "bulge" and "tdisk" for all the quantities 
describing the disk, bulge and thin disk. 

The emissivity of the old stellar population is described 
by an exponential disk and a de Vaucouleurs bulge: 



??(A,J?,2) =7?'^'"^(A,0,0)exp 



R 

■ /jdisk 



ydisk 



B 



+ rj^^'^s^X, 0, 0) exp(-7.67 B^/*) 



^yR^ + {a/by 

Re 



■7/8 



(1) 



(2) 



where R and z are the cylindrical coordinates, ?7^''^'^(A, 0, 0) 
is the stellar emissivity at the centre of the disk, h'^^^^, 
z^'^^ are the scalelength and scaleheight of the disk, 
^buigo(^^ 0, 0) is the stellar emissivity at the centre of the 
bulge. Re is the effective radius of the bulge, and a and b 
are the semi-major and semi-minor axes of the bulge. 

The dust associated with the old stellar population is 
also described by an exponential disk: 



K^-nA,i?,z)=Cf(A,0,0)exp - 



R 



ydisk 



(3) 



where Kgx?(''^; 0, 0) is the extinction coefficient at the cen- 
tre of the disk and h'^'^^ and z^^^^ are the scalelength and 
scaleheight of the dust associated with the old stellar disk. 



The assumption that the scaleheight of the old stellar 
population and associated dust is independent of galacto- 

centric radius is worthy of some comment. For an isother- 
mal disk, the scaleheight, z, is related to the velocity dis- 
persion ^/{v'^), and the density of matter, p, (both at the 
mid-plane) by z oc 

((w2)/p)i/2 (Kellman 1972). One can 
attempt to use this equation to infer the variation of the 
scaleheight of the dust with galactocentric radius by con- 
sidering the gaseous layer in which the dust might be pre- 
sumed to be embedded. Since the velocity dispersion of 
this gas has a value of about lOkms"^, and is observed 
to remain approximately constant with increasing radius 
(van der Kruit & Shostak 1982, Shostak & van der Kruit 
1984, Kim et al. 1999, Sellwood & Balbus 1999, Combes & 
Becquaert 1997), the above equation would imply that in 
an exponential disk in which the matter distribution fol- 
lows that of the stellar luminosity, the scaleheight of the 
gas should slowly increase with radius, with a character- 
istic scalelength twice that of the exponential stellar disk. 
However, when the additional gravitational force of the 
gas and dark matter halo is taken into account (Narayan 
& Jog 2002), the increase of scaleheight with radius is 
predicted to be reduced. Also, the dust-to-gas ratio could 
well be a decreasing function of z, further reducing the 
radial dependence of the dust scaleheight. For the case 
of the stars, it is known that the velocity dispersion falls 
exponentially with galactocentric radius (Bottema 1993), 
reducing any increase of scaleheight with radius. Therefore 
in the current generation of models we have not attempted 
to make refinements to our assumption that the scale- 
heights of stars and dust are independent of galactocentric 
radius. 

For NGC 891 and other edge-on galaxies, Xilouris et 
al. (1999) derived the geometrical parameters in Eqs. 1- 
3 by fitting resolved optical and NIR images with sim- 
ulated images produced from radiative transfer calcula- 
tions. They found scalcheights for the old stellar popula- 
tion of several hundred parsec, a result which had already 
been known for the Milky Way, and which can be physi- 
cally attributed to the increase of the kinetic temperature 
of stellar populations on timcscales of order Gyr due to 
encounters with molecular clouds and/or spiral density 
waves (Wielen 1977). Another result to emerge from this 
work was that the old stellar populations have scaleheights 
larger than those of the associated dust. The resulting pa- 
rameters for NGC 891, derived independently at each op- 
tical/NIR wavelength by Xilouris et al. (1999), were used 
in Paper I for the modelling of the UV-submm SED of 
this galaxy and were found to be consistent with the FIR 
morphology of NGC 891 (Popescu et al. 2004). Therefore 
we adopt these parameters in this work, taking their av- 
erages over the optical/NIR range, since no trend with 
wavelength is apparent. The exception is the scalelength 
of the old stellar disk /i^ which decreases with increas- 
ing wavelength. The averaged values describing the distri- 
bution of the old stellar population and associated dust 
are given in Table 1, where the length parameters have 
been normalised to the value of /tj'^*^ in the B band. The 
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Table 1. The parameters of the model. All length pa- 
rameters are normalised to the B-band scalelength of the 
disk. 



~disk 

disk 
«d 



[^tdisk 
_tdisk 



Re 

b/a 



f.disk 

f.tdisk 
LB 



0.074 
1.406 
0.048 



1.000 
0.016 
/i*/'"'' 1.000 

„tdisk 
«d 



0.016 



0.229 
0.6 



0.387 



Table 2. Wavelength dependence of the scalelength of the 
disk normalised to its value in the B band. 



UV B 


V 


I 


J 


K 


hf^'' - 1.000 


0.966 


0.869 


0.776 


0.683 



^tdisk equated to the scalelength of the old stellar pop- 
ulation in the B band (see Table 1). The dust associated 
with the young stellar population was fixed to have the 
same scalelength and scaleheight as for the young stellar 
disk, namely hf'''^=hl'^'^^ and 4<i'*=4<i'* (see Table 1). 
The reason for this choice is that our thin disk of dust was 
introduced to mimic the diffuse component of dust which 
pervades the spiral arms, and which occupies the same 
volume as that occupied by the young stars. This choice 
is also physically plausible, since the star-formation rate 
is closely connected to the gas surface density in the spi- 
ral arms, and this gas bears the grains which caused the 
obscuration. In principle, we might expect that the metal- 
licity gradient of the gas within the galaxy would decrease 
the radial scalelength of the dust (/ij^'^"''^) below that of the 
gas. However, the ratio of the gas to stellar surface densi- 
ties increases with galactocentric radius, and so tends to 
cancel any variation in the ratio of the radial scalelength 
of the dust with respect to that of the young stars. 



adopted values for /ig'^*^ (again normalised to the value 
of hg^^^ in the B band) are given in Table 2, where the I 
band value has been interpolated from the J and K band 
values. We should stress that the normalised parameter 
values from Table 1 arc close to the average values found 
by Xilouris et al. (1999) for his sample of edge-on galaxies. 

The emissivity of the young stellar population is also 
specified by an exponential disk: 

where 7/*'^"^'^(A, 0, 0) is the stellar emissivity at the centre 
of the thin disk and /i*'^'^'^ and zl'^^^^ are the scalelength 
and scaleheight of the thin disk. 

The dust associated with the young stellar population 
is again specified by an exponential disk: 

KlftiX, R, z) = 4r (A, 0, 0) exp (-^ - (5) 

where Kg^'|^(A, 0, 0) is the extinction coefficient at the cen- 
tre of the thin disk and h^^^^^ and z^J^^^^ are the scalelength 
and scaleheight of the dust associated with the young stel- 
lar disk. To minimise the number of free parameters, we 
fix the ratio nf^H^ 0, ^)/fAxf'{\ 0, 0) to the value found 
for NGC 891 in Paper I, and express it in terms of the ra- 
tio 7-^'disk^^^,tdisk ^£ central face-on optical depths in the 
B band (see Table 1). This is also close to what was found 
for NGC 5907 (Misiriotis et al. 2001). 

The young stellar population is known to have smaller 
scale heights than the older stellar population and asso- 
ciated dust, and to be seen towards more optically thick 
lines of sight. Typically this population cannot be con- 
strained from UV images, and therefore its scalcliciglit 
_j,tdisk ^jj^g gjjnpiy fixed to be 90 pc (close to that of the 
Milky Way; Mihalas & Binney 1981) and its scalelength 



2.2. The clumpy component 

By their very nature, star-forming regions harbour opti- 
cally thick clouds which are the birth places of massive 
stars. There is therefore a certain probability that radia- 
tion from massive stars will be intercepted and absorbed 
by their parent clouds. This process is accounted for by a 
dumpiness factor F which is defined as the total fraction 
of UV light which is locally absorbed in the star-forming 
regions where the stars were born. Astrophysically this 
process arises because at any particular epoch some frac- 
tion of the massive stars have not had time to escape the 
vicinity of their parent molecular clouds. Thus, F is re- 
lated to the ratio between the distance a star travels in 
its lifetime due to it's random velocity and the typical di- 
mensions of star-forming complexes. To conclude, in our 
formulation the clumpy distribution of dust is associated 
with the opaque parent molecular clouds of massive stars, 
and is not to be confused with the randomly distributed 
clumps (of dust unrelated to the stellar sources). 

The attenuation by the clumpy component has a dif- 
ferent behaviour than that of the diffuse component. One 
difference is that the attenuation by the clumps is indepen- 
dent of the inclination of the galaxy. Another difference is 
that the wavelength dependence is not determined by the 
optical properties of the grains (because the clouds are so 
opaque that they block the same proportion of light from 
a given star at a given time at each wavelength), but in- 
stead arises because stars of different masses survive for 
different times, such that lower mass and redder stars can 
escape further from the star-forming complexes in their 
lifetimes. A proper treatment of the dumpiness factor is 
important since the dumpiness will change the shape and 
the inclination dependence of the UV attenuation curves 
of star-forming galaxies. 
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Table 3. The wavelength dependence of the normalised 
extinction coefficient k*^^ (normalised to its value in the 
B band), albedo and the anisotropy factor g. 
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2.44 


0.43 


0.53 


2200 


2.60 


0.43 


0.49 


2500 


1.99 


0.52 


0.48 


2800 


1.66 


0.54 


0.49 


4430 


1.00 


0.53 


0.50 


5640 


0.74 


0.52 


0.45 


8090 


0.44 


0.49 


0.36 


12590 


0.20 


0.37 


0.15 


22000 


0.06 


0.16 


0.04 



^The wavelengths are those for which the attenuation is calcu- 
lated (see Sect. 3). 



2.3. The dust model 

The dust model corresponds to the graphite/silicate mix 

of Laor & Draine (1993) and to the a~^^^ grain size dis- 
tribution of Mathis, Rumple & Nordsieck (1977) and was 
also used in the calculation of the optical-submm emission 
in Paper I. The wavelength dependence of the extinction 
coefficient Kext) albedo and the anisotropy factor g for this 
model is given in Table 3. This model is consistent with 
a Milky Way extinction curve and, as already noted, with 
the extinction curves found for the galaxies modelled in 
Paper II. 

2.4. The free parameters of the model 

The optical-submm SED model presented in Paper I has 
only three free parameters, since, as described before, the 
geometry is constrained by optical/NIR images. The three 
free parameters are: SFR, mass of dust in the thin disk, 
and the dumpiness factor F. For the specific application of 
calculating attenuation in the UV to NIR spectral range, 
SFR is not needed, since extinction does not depend on 
the strength of sources. Also, because here we have fixed 
the ratio between the Tg'^'^'^ and T-g'^'^*^ (see Table 1), and 
because the dust model is also fixed, the mass of dust in 
the thin disk is fully determined b^ the total central face- 
on optical depth in B band, = rg'^^^'^-l-TB'*'^'^^. The third 
free parameter for the calculation of the whole SED - the 
dumpiness factor F - remains a free parameter for the 
calculation of attenuation. In addition, attenuation also 
depends on the inclination i. In summary, we need three 
parameters to fully determine the attenuation in a galaxy: 
Tg, F and i. A fourth free parameter - the bulge-to-disk 
ratio - is introduced in Sect. 5 to account for the different 
morphologies encountered in the Hubble sequence of spiral 
galaxies. 



3. Calculation of attenuation in the diffuse 
component 

This section describes the use of radiation transfer cal- 
culations to obtain the attenuation in the diffuse com- 
ponent. No radiation transfer calculations are needed for 
the clumpy component, which is handled analytically (see 
Sect. 5). 

The basic approach used here is to calculate the atten- 
uation separately for the three diffuse geometrical compo- 
nents of our model: the disk, the thin disk and the bulge. 
In each case the same fixed geometry of the diffuse dust 
is adopted, which is the superposition of the dust in the 
disk and the dust in the thin disk. In other words we de- 
rive the attenuation of the old stellar population in the 
disk as seen through the dust in the disk and thin disk; 
we derive the attenuation of the young stellar population 
in the thin disk as seen though the same dust in the disk 
and thin disk; and we derive the attenuation of the bulge, 
also viewed through the dust in the disk and thin disk. 

The calculations were performed for combinations of 
the two parameters affecting the diffuse component, 
and inclination i. For the sampling in Tg we chose the set 
of values: 0.1, 0.3, 0.5, 1.0, 2.0, 4.0, 8.0 which range from ex- 
tremely optically thin to moderately optically thick cases. ^ 
For the sampling in inclination we chose < cos(i) < 1, 
with A cos(i) = 0.05. Each calculation was performed at a 
different wavelength, covering the whole UV/optical/NIR 
range, such that attenuation can be conveniently calcu- 
lated for both standard optical bands and the bands cov- 
ered by GALEX. Our choice of UV wavelengths also sam- 
ples the 2200 A feature. The calculations in the UV range 
were only performed for the thin disk, since this is the only 
component of stellar emissivity emitting in this spectral 
range. In total we used 12 wavelengths, listed in Table 3. 

Simulated images of the pure disk, thin disk and bulge 
were calculated for each combination of parameters using 
the radiative transfer code of Kylafis & Bahcall (1987), 
which includes anisotropic multiple scattering. In total we 
produced 3234 images with a pixel size (equal to the res- 
olution) of 0.0066 of the B-band scalelength sam- 
pled every 5 and 10 pixels in the inner and outer disk, 
respectively. The high resolution of the simulated images 
matches the resolution of the optical images of NGC 891 
used in the optimisation procedure for the derivation of 
the disk/bulge geometry. This choice enables not only high 
accuracy in the derivation of attenuation, but also means 
that the resulting model images (after suitable interpola- 
tion) can be used as template images for comparison with 
observed images of real galaxies. The simulated images 
for the disks extend out to a radius of 4.63 B-band scale- 
lengths h'^^^^, which is equivalent to 3.31 dust scalelengths 

^ We should remark that, since is the opacity through the 
centre of the galaxy, where most of the dust is concentrated, a 
choice of 1.0 for tb actually represents an optically thin galaxy 
over almost all of its area. Even for Tb = 4, which is close to 
that found for NGC 891, less than half of the total bolometric 
luminosity is absorbed by dust. 
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Fig. 1. Simulated dust-attenuated images of the disk, thin disk and bulge seen at 78 degrees inclination, for a total 
central face-on optical depth in the B band Tg = 4. Only the right hand side of each image is displayed, with the centre 
of the galaxy lying half way along the left hand border of each panel. The dust disks which produce the attenuation 
are inclined such that they are closer to the observer on the lower half of each image. All images have been displayed 
with the same spatial scale and have the same colour table. For each pair of images, pixels with the same intensity are 
displayed with the same colour. Each image was scaled such that the brightest pixel on the corresponding unattenuated 
image has the value 1.0. In this way one can see the effect of increased opacity with decreasing wavelength, by comparing 
the left hand panel with the right hand panel, for each geometrical component. The top panels show the disk in the B 
band (left) and K band (right). The smaller extent of the disk as seen in the K band is due to the smaller scalelength 
of the stellar emissivity in the K band compared to that in the B band. In both images the maximum brightness comes 
from the central regions, though at B band the area of bright emission (white colour) is smaller than that in the K 
band, due to the increased opacity. The middle panels show the thin disk at 1650 A (left) and 2800 A (right). In both 
images the maximum brightness is no longer in the central region, but at an intermediate radius, leaving a central 
"hole" , whose area (red colour) is larger at the shorter UV wavelength. The asymmetry around the horizontal axis is 
caused by the stronger forward scattering which brightens the near side of the thin disk compared to the far side. The 
bottom panels show the bulge in the B (left) and K band (right). In the K band bulge emission from both above and 
below the dust disks is seen, though the emission from below is somewhat attenuated. In the B band almost all the 
emission from below the dust disk is obscured. 

^disk^ The simulated images for the bulge extend out to 
1.45 B-band scalelength h^^^^. For both the disks and the 
bulge the extent of the integration is a numerical limit. 



Examples of simulated images for the disk, thin disk and 
bulge are given in Fig. 1. 

We also produced the corresponding intrinsic images 
of the stellar emissivity (as would be observed in the ab- 
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Fig. 2. Examples of the dependence of attenuation (Am) on inclination (i) for the main geometrical components of 
our model: disk (top left), bulgc(top right), and thin disk (bottom left and right). The examples are plotted for the 
B band for the disk and bulge, and both for the B band and UV 1350A for the thin disk. In each panel we plotted 
(from top to bottom) 7 attenuation curves, corresponding to r^: 8, 4, 2, 1, 0.5, 0.3, and 0.1. The face-on orientation 
corresponds to 1 — cos(i) = 0.0 and the edge-on orientation corresponds to 1 — cos(z) = 1.0. 



sence of dust). The attenuation Am was then obtained by 
subtracting the integrated magnitude of the dust affected 
images from the integrated magnitude of the intrinsic im- 
ages. 

At this point we have obtained values of the atten- 
uation for different combinations of i and Tg, at each 
wavelength, and independently for the disk, thin disk and 
bulge. To facilitate access to this information and to allow 
attenuation to be calculated for any inclination, we fit the 
attenuation curves (Am vs i) with polynomial functions 
of the form: 



J2 aj(l - cos(z))j for 1 - cos(z) < 0.90 



Am = <! i=°'" 
bo 

h 



for 1 
for 1 



cos(i) 
cos(z) 



0.95 
1.00 



(6) 



where n = 5 for the disk and thin disk and n — A for the 
bulge. Values for the coefficients aj and for the constants 
ho and hi are given in Tables 4, 5 and 6 for the disk, thin 
disk and bulge, respectively. To derive the attenuation for 
any desired combination of i and Tg one can first apply 
Eq. 6 to obtain the values of Am corresponding to the 
sampled Tg and then interpolate between these values. 
Subsequently one can also interpolate in wavelength. 

It should be noted that slight negative values for Am 
are obtained for some combinations of low inclination and 
low Tg. This mild amplification is due to the scattering of 
light which removes photons travelling at high inclinations 
(in the plane of the disk) and sends them into directions 
with low inclinations, as also shown by Baes & Dejonghe 
(2001) for the case of isotropic scattering. 
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In Fig. 2 we show examples of the dependence of atten- 
uation on incUnation for the three geometrical diffuse com- 
ponents, both in the optical and in the UV. As expected, 
the increase in attenuation with increasing inclination is 
stronger for larger than for lower r^, irrespective of 
geometry or wavelength. 



4. Attenuation characteristics for the disk, thin 
disk and bulge 

4.1. Attenuation of the disk 

In Fig. 3 (upper four panels) we show the wavelength de- 
pendence of the attenuation of the disk, from the B band 
to the K band. As expected, the overall level of the at- 
tenuation increases with increasing . Another feature is 
the bunching of the curves at low inclinations and low r^, 
followed by an increase in the spacing of the curves when 
proc;ecding to higher inclinations and higher Tg. This can 
be explained as follows. At low inclination and Tg most 
of the disk is optically thin and the attenuation will scale 
as — Zo.g(l — r), where t is the line of sight optical depth. 
Thus, in this regime, changes in inclination or induce 
only small changes in apparent magnitude. However, as we 
proceed to higher inclinations and Tg , the disk area which 
is optically thick will increase from the centre, spreading 
to the outside. In the inner optically thick part of the disk 
attenuation scales linearly with t whereas in the outer 
optically thin parts the attenuation still scales logarithmi- 
cally. Thus the increasing fraction of the disk area which 
is optically thick will have the effect of inducing bigger 
and bigger changes in the apparent magnitudes. 

The bunching of the curves of attenuation versus wave- 
length at low inclination and low Tg can also be appreci- 
ated from examination of the curves of attenuation versus 
inclination shown in Fig. 2a (for the case of the B band). 
The bottom curves of Fig. 2a, corresponding to low Tg, 
are almost flat over most of the inclination range, mean- 
ing that the increase in inclination produces only small 
changes in attenuation. Similarly, the increase in the spac- 
ing of the curves of attenuation versus wavelength at high 
inclination and high Tg from Fig. 3 (upper four panels) is 
reflected by the upper curves of Fig. 2a (corresponding to 
high Tg), which continuously steepen. 

A further feature of the disk attenuation is the smaller 
spacing of the curves between the edge-on geometry and 
the slightly inclined from edge-on geometry (the two top 
curves in each of the upper four panels of Fig. 3), an effect 
which becomes more pronounced at high Tg. This effect 
can also be seen in Fig. 2a for the case of the B band, 
where the curves of attenuation versus inclination flat- 
ten when approaching the edge-on viewing angle. This is 
because the stellar disk has a higher scale height than 
the dust. In the edge-on view, the high z tail of the stel- 
lar population will be visible both above and below the 
plane. This will tend to boost the apparent luminosity of 
the disk compared with the slightly inclined disk, partially 



cancelling the dimming due to the increased line-of-sight 
optical depth. 

4.2. Attenuation of the bulge 

Analogous to the case of the disk, we present the wave- 
length dependence of the attenuation of the bulge in Fig 3 
(lower four panels). The most interesting result is the 
larger overall attenuation of the bulge, as compared to 
that of the disk (Fig 3, upper four panels). This is be- 
cause the bulge is more concentrated towards the opti- 
cally thick part of the dust disks. For the same reason the 
curves also exhibit a strong dependence of attenuation on 
inclination. Especially for low Tg, this dependence is even 
more pronounced than for the disks. At high Tg though, 
the variation of attenuation with inclination is limited by 
the larger vertical extent of the stellar population of the 
bulge compared to the vertical extent of the dust. 

Another feature of the attenuation curves for the 
bulges is the crossing of the two uppermost curves, for 
87 and 90 degrees (Fig. 3, lower four panels). The inter- 
section of the curves occurs at a critical wavelength (which 
depends on Tg), shortwards of which the highest attenua- 
tion is no longer for 90 degrees inclination, but instead at 
a slightly smaller inclination. This critical wavelength in- 
creases with increasing Tg: it lies between B and V bands 
for Tg = 1 , between V and I bands for Tg = 2, between I 
and J bands for Tg = 4 and goes to K band for Tg = 8. The 
explanation for these features is as follows. If we consider 
that the dust disk is optically thick for all lines of sight 
intersecting the bulge, then the fraction of light blanked 
out is simply proportional to the projected area of the disk 
onto the bulge. Since the projected area of the disk is at 
a minimum for an inclination of 90 degrees, moving away 
from the edge-on geometry will increase the attenuation. 
If we progress to even lower inclinations the disk will start 
to become optically thin for some lines of sight through 
the bulge, and then the attenuation will start to decrease 
with decreasing inclination. This non-monotonic progres- 
sion of attenuation with inclination close to the edge-on 
orientation is also illustrated in Fig. 2b (for the exam- 
ple of the B band). There one can see that the maximum 
attenuation is no longer reached at 90 degrees for Tg > 1. 

4.3. Attenuation of the thin disk 

In Fig. 4 (upper four panels) we show the dependence 
of the attenuation of the thin disk on UV wavelength. 
In this wavelength range the opacity does not decrease 
monotonically with increasing wavelength, as it does in 
the optical regime, because of the local maximum around 
2200 A in the extinction efficiencies. This causes the bump 
around 2200 A seen in the attenuation curves in Fig. 4 
(upper four panels). 

As in the case of the attenuation of the disk (Fig. 3, up- 
per four panels), one can observe a bunching of the curves 
for low inclinations and low Tg, followed by an increase in 
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Fig. 3. The dependence of the attenuation (Am) of the disk (upper four panels) and of the bulge (lower four panels) 
on optical/NIR wavelength (A). In each case the panels display the attenuation curves for different Tg: 0.3, 1, 4 and 
8. In each panel the symbols represent attenuations calculated for 5 optical/NIR wavelengths: B, V, I, J and K band. 
The solid lines represent the linear interpolation between adjacent calculated points. In each panel we plotted (from 
top to bottom) 8 attenuation curves, corresponding to inclinations 1 — cos(i): 1, 0.95, 0.90, 0.85, 0.75, 0.60, 0.40, 0. 
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the spacing between curves in proceeding to high inclina- 
tions and high . In addition to the factors contributing 
to this effect discussed in Sect. 4.1, scattering of light also 
payes a role. Photons travelling at high inclinations (in the 
plane of the disk) are scattered into directions with low 
inclinations, such that they can escape the disk (as ex- 
plained by Baes & Dejonghe 2001 for the case of isotropic 
scattering). This factor also influences the attenuation in 
the disk, but is more important in the case of the thin 
disk, since here the stars have a stronger correlation with 
the dust, so that the contrast between the edge-on and 
the face-on optical depth is higher. 

At high Tg, however, there is a small tendency for the 
curves for the edge-on and nearly edge-on inclinations to 
come together. This effect, which is most prominent for 
912 A and = 8 (Fig. 4, upper four panels), is primarily 
a saturation effect due to the fact that only a thin skin is 
seen close to edge-on. 

In Fig. 4 (lower four panels) we show the dependence 
of the attenuation of the thin disk on optical wavelength. 
At larger inclinations, the thin disk exhibits the largest 
attenuation in the optical range, as compared with the 
disk and the bulge. This is because of the strong spatial 
correlation between the stars and the dust. For the same 
reason, the thin disk also exhibits the strongest depen- 
dence on inclination. Despite the larger attenuation, the 
contribution of this geometrical component to the overall 
attenuation of the galaxy in the optical range is small and 
can be neglected, since the thin disk emits mainly in the 
UV. However the attenuation of the thin disk in the op- 
tical range is important for the calculation of optical line 
emission arising from the thin disk. 

4.4. Intercomparison of the attenuations of the disk, 
thin disk and bulge 

In the previous subsections we described and discussed the 
main characteristics of the attenuation curves of the main 
geometrical components of our model, namely of the disk, 
thin disk and the bulge. To make an intercomparison of the 
curves for the different geometries we superpose in Fig. 5 
attenuation curves for the three geometrical components, 
for the same inclination and Tg . The value of Tg was chosen 
to be 4, which, for the geometry of our model, is close to 
what we found for NGC 891. 

At low and intermediate inclinations the bulge suffers 
a larger attenuation than that of the disk or thin disk 
(Fig. 5a, b). As already noted, this is because the bulge is 
more concentrated towards the optically thick part of the 
dust disks. In addition, and for the same reason, the figure 
also shows that the attenuation of the bulge has the steep- 
est dependence on optical wavelength. It is also apparent 
from Fig. 5a,b that the attenuation curves for the disk 
and thin disk show quite similar behaviour to each other, 
indicating that at low and intermediate inclinations the 
curves are not sensitive to the differences in scale height 
and scale length between the two disks. There is neverthe- 
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less a small difference between the curves for the disk and 
thin disk, in the sense that the thin disk suffers slightly 
less attenuation, except for the B band. This is because 
the disk has smaller scale lengths for the stellar popula- 
tion than the thin disk, except for the B band, where the 
scale length is the same. 

With increasing line of sight optical depth, which is 
the same as increasing inclination and/or decreasing wave- 
length, the attenuation curve for the thin disk starts to di- 
verge from that of the disk (Fig 5c); it will increase until 
it intersects with the attenuation curve for the bulge. This 
is because the thin disk has the smallest scale height, the 
strongest spatial correlation with the dust and the largest 
contrast between the face-on and edge-on orientation, and 
so its attenuation depends most strongly on inclination 
(or wavelength). So for edge-on orientation (Fig. 5d) the 
attenuation of the thin disk dominates that of the other 
geometrical components, reversing the situation for the 
face-on orientation (Fig. 5a), where the thin disk has the 
smallest attenuation. 

5. Composite attenuation curves for spiral 
galaxies 

Spiral galaxies have varying proportions of their stellar 
luminosities emitted by the disk, thin disk and bulge 
and also varying amounts of dumpiness. For instance the 
bulge-to-disk ratio decreases in going from earlier to later 
spiral types (see for example Fig. 7 of Trujillo et al. 2002). 
In the previous sections we have calculated attenuation 
curves quantifying the extinction of the individual stellar 
populations from the disk, thin disk and bulge due to the 
diffuse dust. These can be combined to obtain a compos- 
ite attenuation curve for the light illuminating the diffuse 
dust in a galaxy. If we also consider the photons that don't 
reach the diffuse dust because they arc absorbed by the 
clumpy dust, local to the star-forming regions, then we 
obtain the composite attenuation curve for the total lumi- 
nosity of the galaxy. 

5.1. A formula to calculate composite attenuations 

At a given wavelength A, the total attenuation Amx in a 
galaxy is given by: 

Amx = -2.5 log ^, (7) 

where 5° and Sx are the intrinsic and the apparent flux 
densities, respectively. The quantities and Sx can be 
further expressed as a summation of the corresponding 
quantities for the disk, thin disk and bulge: 

^0 ^Ojdisk . ^Ojtdisk , ^0, bulge /g\ 

Sx = Sf'^ + Sl^''^ + Sl''^^\ (9) 

The apparent and intrinsic flux densities for the disk, thin 
disk and bulge are related as follows: 

^o,disk ^ 5disk;^Q 2.5 (10) 
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Fig. 5. Comparison between the attenuation curves of the disk, thin disk and bulge. The 4 panels show the dependence 
of the attenuation on optical wavelength, for 4 inclinations: 0, 60, 87 and 90 degrees, and for the same = 4. In each 
panel the attenuation for the disk, thin disk and bulge are plotted with solid, dashed and dotted lines, respectively. 



r*0,tdisk 



cftdisk 

10 



Ffx 



2.5 



^O.bulge ^ S'bulge-^Q ^ 



(11) 
(12) 



where Am^"''* 



, Am 



tdisk Am'?"'^'' are the attenuation 



values for the disk, thin disk and bulge, which can be de- 
rived from Tables 4-6, F is the dumpiness factor as defined 
in Sect. 2.2 and f\ is the function defined in Appendix A 
and tabulated in Table A.l. In physical terms, F fx is the 
fraction of the emitted UV flux density at wavelength A 
which is locally absorbed in star-forming regions. 

If we define the ratios of the apparent luminosities of 
disk, thin disk and bulge to the total apparent luminosity 
to be: 



^disk 



^disk 



(13) 



Ctdisk 
tdisk _ "^A 

^ Sx 

rrbulgc 


Sx ■ 



bulge 



(14) 
(15) 



then the total attenuation of the galaxy can be expressed 
as: 



Amx = 2.51og(r^'"''10 



^disk 



bulge 



+ 



Ffx 



Am\^'''^ 
10 2.5 



Am^"'*^'' 
2.5 



(16) 
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This is a general formula giving the composite attenua- 
tion curve as a function of 5 independent parameters: i, 

and F. For practical apphcation, this for- 



-f „disk -,bulgo 



' A 



A 



mula can be simplified if we write it separately for the UV 
and the optical/NIR ranges. It is a good approximation 
to consider that the thin disk emits only in the UV range 
and the disk and bulge only in the optical/NIR range. This 
approximation has been validated by our detailed analysis 
on NGC 891 in Paper I. In the UV Eq. 16 becomes: 



2.5log{l-Ffx), 
while in the optical/NIR Eq. 16 becomes: 



^^optical 



2.5l0g((l 



bulge 

'■a 



^disk 



) 10^:3- 



+ 



^bulge 



2.5 ). 



(17) 



(18) 



method to derive Tg would be to use the emission from 
the Balmer recombination lines, integrated over the whole 
galaxy. In our model the emission from these lines arises 
from the thin stellar disk. These lines will be attenuated 
both by the diffuse dust (associated with both the young 
and old stellar populations) and by the clumpy dust (as- 
sociated with the star-forming regions). Since in the for- 
mulation of our model the fraction of the emission locally 
absorbed in the clumpy component is the same for all lines 
in the Balmer series, the line ratio will depend only on the 
attenuation by the diffuse dust. This attenuation can read- 
ily be found by interpolating in wavelength between the 
attenuation values for the thin disk derived from Table 5. 
As an example, the ratio Ra/H/S was derived by interpo- 
lating in wavelength between the adjacent optical bands to 
the wavelength of the recombination lines. To present the 
line ratio in an analogous way to the information from 
Tables 4-6, we fitted the log(Ha/H/3) with a 5th order 
polynomial function of the form: 



5.2. Choice of parameter values 

Both Eqs. 17 and 18 have three free parameters: r^, i 
and F for the UV range and Tg, i and r^"'^° for the opti- 
cal/NIR range. 

For resolved objects, the apparent bulge-to-total lumi- 
nosity ratio r^"'^° can be derived directly from observa- 
tions, by decomposing the bidge from the disk. 

A realistic derivation of the F factor requires a com- 
plete modelling of the whole UV-submm SED. Another 
approach would be to derive the F factor directly from 
comparing highly resolved FIR images with images in the 
UV, which soon may become feasible with the advent of 
GALEX in conjunction with ISO and future SIRTF obser- 
vations of nearby galaxies. For a simple recipe, the median 
value of the F factors derived in Paper II can be used, 
namely F — 0.22, which is also the value obtained for 
NGC 891. We note here that it cannot be ruled out that F 
may increase with SFR, as suggested by the study of the 
FIR/radio correlation by Pierini et al. (2003). Physically, 
this might be expected if an increased SFR is accompanied 
not only by an increase in the number of independent HII 
regions, but also by a higher probability for further star 
formation to happen preferentially near already existing 
HII regions. Then, the F factor would also increase, as a 
consequence of the increased blocking capability of the op- 
tically thick molecular clouds in the star-forming complex. 
This would be expected to occur if star formation is a self 
propagating phenomenon, in which preceding generations 
of stars can trigger the formation of new generations. 

The values of opacities Tg are also best found from a 
complete modelling of the whole UV-subrnni SED. In a 
companion paper we will present a grid of calculations for 
the FIR/submm emission corresponding to the attenua- 
tions presented in Tables 4-6, from which Tg can be ex- 
tracted from a self-consistent fit to the entire UV-submm 
range. In the absence of FIR photometry, an alternative 




1 - cos(i) < 0.90 

1 - cos(i) = 0.95 
1 - cos(i) = 1.00 



(19) 



where n = 4. Values for the coefficients aj and for the con- 
stants 60 and bi are given in Table 7 and the dependence 
of the Hq;/H/3 line ratio on inclination and Tg is shown 
in Fig. 6. If a measurement of the line ratio Ra/ILp were 
available for a galaxy of known inclination, this could be 
compared with the calculated ratios derived from Table 7 
to derive Tg. 



Table 7. Coefficients and constants for the calculation of 
the Ha /up line ratio 
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-0.042 


2.435 


-4.730 


3.709 
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5.3. An example 

To illustrate solutions derived from Eqs. 17-18, we calcu- 
late composite attenuation curves for NGC 891 as would 
be obtained for different inclinations. For this galaxy 
Tg = 4.2 and F = 0.22 (see Paper I). In this particular 
case the apparent bulge-to-total ratio r^"*'' can be de- 
rived from a known intrinsic value (0.246, averaged over 
the optical/NIR; Xilouris et al. 1999) by application of 
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Fig. 6. The effect of dust on the Ha/H/3 line ratio as a 
function of inclination i. The seven curves (from top to 
bottom) correspond to r^: 8, 4, 2, 1, 0.5, 0.3, and 0.1. The 
face-on orientation corresponds to 1 — cos(i) = 0.0 and the 
edge-on orientation corresponds to 1 — cos(i) = 1.0. 
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Fig. 7. Composite attenuation curves for the galaxy 
NGC 891, as would be obtained if the galaxy were viewed 
at inclinations of 0, 60 and 90 degrees, = 4.2 and 
F = 0.22. 



Eqns. 10 and 12 for the wavelengths and inclinations of 
interest. Using the values tabulated in Tables 4-6 we ob- 
tained the curves plotted in Fig. 7 for three inclinations: 
0, 60 and 90 degrees. The curves show a smooth progres- 
sion with wavelength, except for the edge-on orientation, 
where there is a steep step between optical and UV, due 
to the large contrast in attenuation between the thin and 
thick disks viewed edge-on. 

Using the attenuation curve for NGC 891 {i = 89.8°) 
and the intrinsic flux densities from Paper I we can pre- 
dict the apparent magnitude of the galaxy in the UV and 
compare this prediction with observations. Until now the 
galaxy has been detected only at 2500 A by Marcum et 
al. (2001) who measured m^^oo = 13.00 ± 0.23. Our pre- 
diction is m^"'''^*' = 13.29, in good agreement with the ob- 
servations, thus reinforcing the validity of the model for 
the intrinsic stellar luminosity and attenuation presented 
in Paper I. 

6. The effect of the bulge-to-disk ratio on the 
derived opacities 

To illustrate the use of composite attenuation curves in 
the interpretation of optical/NIR data we consider the 
effect of a varying bulge-to-disk ratio on the inclination 
dependence of the apparent optical emission from galaxies. 
Traditionally, astronomers have used this dependence to 
make statistical estimates of the opacity of galaxian disks 
through analysis of large samples of spiral galaxies. But 
these studies have not taken into account the fact that 



the attenuation of the bulge has a different variation with 
inclination than that of the disk, as shown in Sect. 4. 

In order to quantify the effect of the bulge-to-disk ra- 
tio on the derived opacities, we calculated the variation of 
the composite attenuation with inclination for three dif- 
ferent opacities, = 2,4,8. These curves were calculated 
for a grid of values of the bulge-to-disk ratio ranging from 
to 100, which represent the whole range between "pure 
disk" and "pure bulge" galaxies. We chose to do these 
calculations in the K band in order to compare our re- 
sults with recent statistical studies of internal extinction 
in spiral galaxies in the NIR. For a convenient compari- 
son with observational studies we fitted the variation of 
the composite attenuation with -log(cos(i)) using a linear 
fit^ to the lower inclination range (1 — cos(i) < 0.7). The 
resulting slope of the fit 7k is plotted in Fig. 8 versus the 
bulge-to-disk ratio and for the 3 opacities considered. All 
curves show a "S-like" shape, tending towards the asymp- 
totes of the "pure" disk and "pure" bulge. It can be seen 
from Fig. 8 that an increase from to 1 in bulge-to-disk 
ratio (which embraces the observed range; Trujillo et al. 
2002) produces a comparable change in 7k to that result- 
ing from a doubling of the opacity of a "pure" disk. For 
example, changing the bulge-to-disk ratio from to 1 for 
Tg = 2 increases 7k by a factor of 2.2, whereas increasing 
the opacity from Tg = 2 to 4 changes 7k by a factor of only 
1.9 for a "pure" disk galaxy. Thus, increasing the bulge- 
to-disk ratio at a constant opacity can mimic the effect of 
increasing the opacity of a "pure" disk. 

^ -log(cos(i)) is equal to log(a/6), where 0/6 is the axial ratio 
of the galaxy. 
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Fig. 8. The slope 7k versus the intrinsic bulge-to-disk ra- 
tio, where 7k is the slope of a linear fit to the low inclina- 
tion part (1 — cos(i) < 0.7) of the composite attenuation 
curve (Am vs -log(cos(i))) in the K band. The calculations 
were performed for Xg = 8, 4, 2. 

This may have consequences for the use of optical sta- 
tistical samples to evaluate the dependence of opacity on 
galaxy luminosity. For example, in their study of 15224 
spiral galaxies from the 2 Micron All-Sky Survey, Masters 
et al. (2003) found a trend for the slope 7k (and also 
for the corresponding slopes in J and H bands) to in- 
crease with increasing K-band luminosity. These authors 
interpret this result as a trend of increasing disk opac- 
ity with increasing K-band luminosity. However, galaxies 
that have brighter K-band luminosities are also biased to- 
wards earlier- type spirals (e.g. Bosclli et al. 1997), which 
in turn have larger bulge-to-disk ratios. Fig. 13 of Masters 
et al. (2003) show a variation of 0.2 in 7k over the lumi- 
nosity range of their sample, which is comparable to the 
variation in 7k found in our Fig. 8 between bulge-to-disk 
ratios of and 1. This raises the possibility that some of 
the trend found by Masters ct al. (2003) is due to a sys- 
tematic increase in the bulge-to-disk ratio with luminosity 
in their sample. It will be important to quantify this effect 
in studies that investigate the dependence of the ratio of 
the FIR luminosity to the intrinsic UV luminosity of gas 
rich galaxies on the stellar mass of the galaxies (Pierini & 
MoUer 2003). In general, ignoring the presence of bulges 
can lead to a systematic overestimate of the opacity of 
disks. 

7. Summary and conclusions 

Wc present new calculations for the attenuation of the 
integrated stellar light from spiral galaxies, utilising ge- 
ometries for stars and dust constrained by a joint con- 



sideration of the UV/optical and FIR/submm SEDs. In 
addition to the single exponential diffuse dust disk used 
in previous studies of attenuation, wc also invoke a clumpy 
and strongly heated distribution of grains spatially corre- 
lated with the star-forming regions (required to account 
for the FIR colours) and diffuse dust associated with the 
spiral arms (needed to account for the amplitude of the 
submm emission). The latter is approximated by a second 
exponential dust disk having the same spatial distribution 
as the young, UV emitting stellar population, also approx- 
imated by an exponential disk - the "thin disk" . The old, 
optical/NIR-emitting stellar population is specified by an 
exponential disk of larger scale height - the "disk" , plus a 
de Vaucouleurs bulge - the "bulge" . 

Radiation transfer calculations were performed sepa- 
rately for the three diffuse emissivity components (disk, 
thin disk and bulge) seen through the same fixed distribu- 
tion of diffuse dust. Wc used the radiative transfer code of 
Kylafis & Bahcall (1987) which includes anisotropic multi- 
ple scattering. The attenuation of each component was cal- 
culated for a grid of central face-on B-band optical depth 
T^'^^^^ and inclination i, for a range of wavelengths from 
912 A to 2.2 yum. The resulting curves of attenuation ver- 
sus inclination were fitted with polynomial function and 
the coefficients tabulated in Tables 4-6. The tables allow 
the attenuation to be obtained for any desired combina- 
tion of Tg^'^*^, i, and wavelength, by interpolation in t^'^^^^ 
and wavelength. 

The clumpy component was handled analytically, 
whereby the wavelength dependence of local absorption 
of starlight in star-forming regions was purely determined 
from geometrical considerations. By combining the global 
attenuation of the individual stellar populations from the 
disk, thin disk and bulge due to the diffuse dust with 
the local attenuation due to the clumpy dust associated 
with the star-forming regions, wc obtain a general for- 
mula for the calculation of composite attenuation of the 
integrated emission from spiral galaxies. As an example of 
this formula, we calculated composite attenuation curves 
for NGC 891 as would be seen at different inclinations. 
Using the attenuation curve for NGC 891 at the actual 
inclination i = 89.8° we predicted the apparent magni- 
tude of this galaxy in the UV and found this prediction 
to be in good agreement with the recent observations of 
Marcum et al. (2001). 

We also used our model to derive the ratio Ha/H/3 
as a function of inclination and Tg. This information is 
presented in Table 7 in form of polynomial fits and can 
be used to derive the value of Tg for objects with no 
FIR/submm photometry. 

The detailed analysis of the attenuation properties of 
the individual geometrical components led to the following 
conclusions: 

— For a typical galaxy and in the optical/NIR spectral 
range, the relative attenuation between the disk, thin 
disk and bulge depends strongly on inclination. At low 
and intermediate inclinations, the bulge suffers a larger 
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attenuation than that of the disk or thin disk, and has 
the steepest dependence on wavelength. Furthermore, 
in this regime, the curves of attenuation versus wave- 
length for the disk and thin disk show quite similar be- 
haviour to each other, indicating that the attenuation 
is not sensitive to the differences in scale height and 
scale length between the two disks. At higher inclina- 
tions, however, the attenuation curves for the thin disk 
start to diverge from those of the disk. For the edge-on 
orientation the attenuation of the thin disk dominates 
that of the other geometrical components, reversing 
the situation for the face-on orientation, where the thin 
disk has the smallest attenuation. 
— Increasing the bulge-to-disk ratio at a constant opac- 
ity can mimic the effect of increasing the opacity of a 
"pure" disk. In general, ignoring the presence of bulges 
can lead to a systematic overestimate of the opacity of 
disks. 
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Appendix A: Wavelength dependence of local 
absorption of starlight in star-forming regions 

As introduced by Popescu et al. (2000), the dumpiness factor 
F denotes the total firaction of UV light which is locally ab- 
sorbed in the star-forming regions where the stars were born. 
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Table 4. Coefficients and constants for tlie calculation of 
tlie attenuation of the disk 
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because stars of different masses survive for different times, 
sucli that lower mass and redder stars can escape further from 
the star-forming regions in their Hfetimes. 

We approximate the solid angle O of a parent cloud sub- 
tended at an offspring star of age t with: 

for 
for 



n = 



In the current work we need to understand the wavelength de- 
pendence of the probability for local absorption of UV photons, 
which we denote by Fx , such that 



(A.l) 



where L\ is the intrinsic luminosity density of the galaxy and 
Al and A2 define the UV spectral range. In our formulation Fx 
is determined only by geometric considerations. For a given 
star, the probability of local absorption of a photon is simply 
proportional to the sohd angle of the parent cloud subtended 
at the star. The cloud is considered to be so opaque that a 
photon is absorbed independent of its wavelength if it passes 
through the cloud material. The wavelength dependence arises 



47rj5o 
27rpo[l 



(1 



f ^local ' 
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i ^ ^locali 
t ^ ^local 
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where po is the "porosity factor" giving the typical fraction of 
those lines of sight passing through the cloud which are blocked 
by the cloud material and tiocai is the typical time taken for the 
star to escape from the star-forming region. This simply ap- 
proximates the star-forming regions with porous spheres. The 
porosity of the spheres, averaged over the whole galaxy, is taken 
to be the same for all stellar masses, so that po is independent 
of M. A further approximation is that each star is taken to 
be born at the centre of the cloud and move outwards at a 
constant velocity. 

The expectation value for the fraction of stellar light 
blocked by star-forming regions for stars of zero-age main se- 
quence (ZAMS) mass M, averaged over the lifetime tgai of the 
galaxy is: 



p(M) 



*(tgal - t) $(M) dt 



*o(M) 



in 



/■*gal 

/ 



(A.3) 



*(tgal - t) *(M) dt 



where $ is the initial mass function (IMF), 
formation history of the galaxy and 



to(M) = 



igal -tstar(M) for tstar(M) < tg^l 
for 4tar(M) > tgal 



^ is the star- 



(A.4) 



where istar(A/') is the lifetime of a star of ZAMS mass M. Here 
we consider a Salpeter IMF and an exponentially declining 

star-formation rate: 



= SFR X e 



t/tc 



(A.5) 



where SFR is the current star-formation rate and tc is the time 
constant. 

The power density Px(M) from a star with ZAMS mass M 

is given by: 

/•igal 

Pa(M)= / SEDx{M,t)'^{t^^i-t)^{M)dt (A.6) 

Jto(M) 

Here we approximate the spectral energy distribution 
SED\{M,t) from an individual star with a time independent 
black-body function: 

Bx{T(M)) 



SEDx{M) 



J Bx{T{M))dX 



L(M) 



(A.7) 



where L{M) and T(M) represent the ZAMS luminosity and 
temperature of a star of mass M. 

We derived L{M) from Maeder (1987) and Smith (1983): 
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T(M) was derived from data in Allen (2000), Lang(1980) and 
Schaerer et al. (1996): 



logT(M)= 5^ai[log(^)] 



(A.9) 



where a, = [3.756, 0.601, 0.250, -0.147, -0.078, 0.038]. The de- 
pendence of stellar main sequence lifetime on ZAMS mass M 
was derived from Maeder (1987), Sandage (1957) and Smith 
(1983): 

2.i9(^)-o.430 for ^> 40 
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Fig. A.l. The wavelength dependence of the function f\. 
The solid line indicates the solution for an exponential 
declining star- formation rate with tc = 5 Gyr. The dashed 
line indicates the solution for a constant star-formation 
rate tc = ^c. Both curves were calculated for fgai = 10 Gyr 
and fiocai =3x10'^ yr. 



Thus, if we define a function fx by 
fx = Fxipo = 1) X 



JxjFxipu = l)Lxd\ 
the dumpiness factor F can be expressed as: 



F = 



Fx 
fx 



(A.14) 



(A.15) 



It now remains to calculate /a, which is a function of only 
three parameters: tiocai, tc and tgai. For a typical velocity of 
a star of 5km/s and dimension of star- formation complex of 
150 pc, iiocai is about 3 x 10^ yr, which is adopted here. For 
tgai and tc we take 10 and 5 Gyr, respectively. The resulting 
solution fx is tabulated in Table. A.l, as a function of wave- 
length and is plotted as a solid line in Fig. A.l. The porosity 
factor corresponding to -F = 0.22 is po = 0.368. The assump- 
tions about the star-formation history contained in tc do not 
play an important role in shaping the fx- To check this we also 
calculated fx for the same parameters but tc = oo . This is 
plotted as the dotted line in Fig. A.l. Indeed, there seems to 
be little difference between the two curves. This is to be ex- 
pected, bearing in mind that the timescale for a star to move 
out of the star-forming region is much less than the age of the 
galaxy. In this paper we adopt the solution for tc = 5 Gyr, 
which is the value used for NGC 891 in Popescu et al. (2000). 



Fx is then given by the luminosity weighted average of 
p{M) over the stellar mass distribution: 

Fx = /^;(^)^(^)'^^ (A.11) 

For normal galaxies Fx is completely specified by the porosity 
factor Po, iiocai, tc and tgai. Since 

Fa =poFa(p(, = 1) (A.12) 
Eq. A.l can be rewritten in terms of po as 

^^P^J^F^^^^^ 
C^xdX 
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Table 5. CoefRcients and constants for the calculation of 
the attenuation of the thin disk; Continued 



Table 5. Coefficients and constants for the calculation of 
the attenuation of the thin disk 
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Table 6. Coefficients and constants for the calculation of 
the attenuation of the bulge 
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Table A.l. The wavelength dependence of the function 
fx defined by Eq. A. 13. 
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